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OPTIMUM  THEORETICAL  STRUCTURES  OF  SONAR  SYSTEMS 
EMPLOYING  SPATIALLY-DISTRIBUTED  RECEIVING  ELEMENTS 


By 


F.  Bryn 


ABSTRACT 


The  paper  reviews  the  structure  of  optimum  active  sonar  systems 
employing  spatially  distributed  receiving  elements.  The  systems 
considered  operate  in  the  presence  of  a  time-variable  noise  field 
and  in  a  randomly-varying  transmission  medium.  The  required 
signal-processing  operations  can  be  divided  into  two  groups,  one 
depending  upon  the  noise  characteristics  and  the  other  upon  the 
signal  characteristics.  The  noise-dependent  processor  is  studied 
with  special  emphasis  on  structures  that  can  adapt  themselves  to  a 
time-varying  noise  field.  The  signal-dependent  processor  is 
derived  for  some  simple  types  of  random  signal  distortion.  The 
signal  characteristics  required  for  the  instrumentation  of  such 
processors  are  related  to  the  impulse  response  of  the  target-medium 
combination » 

* 
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1 .  INTRODUCTION 


The  task  of  an  active  or  passive  scnar  system  is  generally  to  survey  a 
finite  ocean  area  in  tho  sense  of  detecting  and  locating  submarines 
moving  into  the  area.  The  task  involves  two  ooerations:  the  acquisition 
of  target  information  from  the  acoustic  field,  and  the  making  of 
decisions  about  the  target  situation  based  upon  the  acquired  information. 
See  Fig.  1.  In  this  paper  we  shall  analyse  the  operation  of  an  active 
system.  The  passive  system  differs  mainly  in  the  typo  of  signals  to  be 
detected.  The  sonar  receiver  is  shown  in  Fig.  2.  It  consists  of  a 
sensing  device  usually  referred  to  as  the  hydrophone  array,  an  information 
processing  device,  and  a  decision  device.  With  the  sensing  device  we 
carry  out  measurements  on  the  random  acoustic  field  in  the  ocean  for  the 
purpose  cf  detecting  acoustic  echo  signals  reflected  from  targets 
within  the  surveyed  ocean  area.  The  information  processing  device  reduces 
the  measurements  to  a  form  required  by  the  decision  process. 

In  this  paper  we  will  describe  the  optimum  structure  of  active  sonar 
systems,  V/e  begin  with  a  study  of  the  ensemble  of  received  echo  signals 
and  indicate  how  this  is  related  to  the  randomness  of  the  transmission 
medium  and  the  targets.  Vie  proceed  with  a  survey  of  the  elements  of 
decision  theory,  and  ocint  out  how  the  optimum  structure  of  the  processing 
device  can  be  deduced  From  the  characteristics  of  the  input  measurements. 

In  many  practical  situations  ws  find  that  the  description  of  the  optimum 
system  depends  on  parameters  which  exhibit  random  variations  in  time. 
Detection  systems  which  measure  or  estimate  the  true  values  of  one  or 
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Surveyed  Ocean 
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Fig.  2  THE  SONAR  RECEIVER 


more  random  parameters  and  adjust  the  structure  of  the  processor  in 
accordance  with  the  estimates,  are  referred  to  as  adaptive  systems.  A 
considerable  part  cf  this  paper  is  allocated  to  a  review  of  techniques 
of  adaption  tc  time-variable  gaussian  ambient  noise  fields. 

The  signal  dependent  part  of  the  signal  processor  is  studied  in 
Chapter  5.  As  is  well  known  the  specific  structure  of  the  processor 
can  only  be  found  for  some  few  types  of  received  signal  characteristics. 
We  study  this  structure  for  the  small  signal  situation  and  for  the  case 
of  signals  composed  of  gaussian  amplitude  perturbations  on  some  mean 
signal.  The  signal  information  required  to  specify  these  processors  is 
discussed  and  related  to  the  target-medium  filters  introduced  in 
Chapter  2. 


2.  SONAR  SIGNALS 

In  a  situation  involving  surveillance  of  submarine  motion  within  an 
ocean  area  by  a  sonar  system,  the  ultimate  information  required  from  the 
system  must  relate  to  the  movements  cf  targets  within  the  surveyed  area*. 
The  most  important  information  concerning  each  target  is  its  position 
and  velocity  relative  to  the  sonar  system  as  functions  of  time.  Relative 
position  is  described  in  terms  of  range  R,  bearing  y  and  depth  d 
whereas  relative  velocity  is  described  in  terms  of  relative  radial 
velocity  v  and  target  aspect  angle  cp  .  Target  parameters  are  referred 
to  a  point  which  is  fixed  relative  tc  the  sonar  system,  e.g.  one- of  the 
hydrophones  in  the  receiving  array.  The  parameters  -will  generally  take 
on  different  values  as  we  move  over  the  aperture  of  the  receiving  array. 
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•It  is  in  fact  this  variation  which  enables  us  to  measure  target  bearing 
and  aspect  angle. 


With  an  active  sonar  system  the  desired  target  information  can  be 
estimated  from  the  temporal  and  spatial  structure  of  the  received  echo 
signals.  Thus  range  is  determined  from  the  absolute  delay  between  a  sonar 
transmission  and  the  received  echo  signal  whereas  bearing  is  determined 
from  the  relative  arrival  times  of  echoes  on  the  hydrophones  of  the 
receiving  array.  Relative  radial  velocity  is  determined  by  the  change 
of  timescale  of  the  echo  signals  from  that  of  the  transmitted  signal  and 
aspect  angle  is  determined  from  combined  measurements  of  range  and 
bearing  of  the  individual  scatterers  of  the  target.  (A  system  with  high 
bearing  and  range  resolution  capabilities  is  required  for  such 
measurements. )  Finally,  target  depth  may  be  estimated  from  the  vertical 
direction  of  incidence  of  the  echo  signal  or  in  some  situations  frem  uhe 
multipath  structure  of  the  echo  signals. 

The  number  of  parameters  that  influence  the  temporal  and  spatial  structure 
of  echo  signals  are  much  greater  than  the  five  mentioned  above.  Soma  of 
these  are  associated  with  the  targets  e.g.  its  orientation,  dimensions  and 
reflecting  properties,  and  some  are  associated  with  the  acoustic  trans¬ 
mission  properties  of  the  ocean  e.g.  the  refracting  and  scattering 
properties  of  the  ocean  volume  and  the  reflecting  properties  of  the  ocean 
boundaries.  We  know  also  that  some  of  these  parameters  are  non-random 
while  others  are  random  timevariable  quantities.  The  group  of  target 
parameters  which  we  wish  to  determine  by  our  sonar  system  will  be 
considered  to  be  the  elements  of  a  target  parameter  vector  a  •  The 
remaining  parameters  which  influence  the  echo  signals  will  be  considered 
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to  ba  elements  of  a  random  vector  r  .  Although  all  the  elements  of  r 
may  not  be  truly  random  in  nature  they  will  be  considered  to  be  so  from 
the  point  of  view  of  the  sonar  receiver  which  must  ba  optimum  on  the 
average  against  these  unknown  parameters.  The  upper  diagram  of  fig.  5 
illustrates  the  important  fact  that  there  are  several  possible  echo 
signals  s  (3.)  for  each  value  of  ql.  In  the  remaining  part  of  this 
chapter  we  shall  introduce  the  concept  of  target  impulse  responses  and 
show  how  these  can  be  used  to  relate  the  target  parameter  vector  a 
tr  a  family  of  received  echo  signals  "s  (a). 


Let  there  be  A  hydrophenes  in  our  receiving  array  and  1st 

h.  (tj,  t^f  a),  i  -  1,  2 .  A,  be  the  echo  signal  received  at 

time  t  on  the  i'th  hydrophone  from  a  target  with  parameters  a  in 

response  to  an  imoulse  emitted  from  the  transmitter  at  time  t^  (fief.  1  ) 

'.Ve  note  that  a  expresses  the  target  parameters  at  the  moment  the  impuls 

impinges  on  the  target.  'Then  measuring  h  (t  ,  t  ,  a  )  we  introduce 

i  1  2 

t  »  +  t  and  study  the  variations  with  t  for  fixed  values  of  t^. 

However,  when  used  in  convolution  integrals  to  find  the  signal  received 
on  the  i'th  hydrcchone  at  time  t^  in  response  to  a  transmitted 
signal  s(t]  we  introduce  t^  =  t^  -  T  and  consider  it  as  a  function 
cf  t  for  fixed  values  of  t  .  This  latter  attitude  will  be  maintained 


throughout  this  chapter.  For  convenience  we  shall  write  t^  3  t, 
t1  a  t  -  t  ,  and  h.  (t  -  r  ,  t,  a  )  3  h^  (  t ,  t,  a).  As  a  function 
of  T  we  assume  h^  (  t  ,  t,  a  )  to  be  the  echo  signal  as  seen  through 
the  input  bandwidth  9  of  the  sonar  system.  As  a  function  of  time  _  t, 
it  is  assumed  to  be  a  stationary  random  function  of  bandwidth  «B.  It 
reflects  the  randomness  residing  in  the  target  and  the  transmission 
medium.  The  set  of  impulse  responses  received  on  the  A  hydrcchones 
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.will  bs  represented  by  the  vector  valued  function 


a  )  =  i  h,  (  t  ,  t,  a  )  j,  i  =  1,  2...  A.  (Eq.  l) 


The  echo  signals  received  from  the  target  in  response  to  a  transmitted 
sonar  signal  s(t)  can  now  be  found  by  convolving  s(t)  with  the 
individual  impulse  responses.  The  set  of  received  signals  will  be 
expressed  by  the  vector  valued  function 


s  (t,  )  = 


1,  2...  A,  (Eq.2) 


where  we  have  .vritten  h  for  h  (  T  ,  t,  a  )»  h^  for  h^  (  <r 
where  s_^  (t,  h.  ;  is  obtained  by  convolving  s(t)  with  h^  (  t 

5i  hi^  =  s(t)  *  lv  (  T  ,  t,  a  ). 


(Eq.  3) 


A  simple  diagram  illustrating  the  formation  of  s^  (t,  h^)  from  s(t] 
is  shown  in  fig.  3a.  For  subsequent  discussions  we  shall  find  it 
convenient  to  change  each  h  (  T  ,  t,  a)  into  two  parallel  connected 

— ft  ml 

filters  g^  (t,  t,  a  )  and  r.  (t  ,  t,  a  ]  as  shown  in  fig.  3fc.  The 
set  of  functions  expressed  by 


t, 


i  =  1 


A,  (Eo.4) 


describe  the  average  impulse  responses  of  the  target-medium  combination 
as  seen  on  the  various  hydrophones.  Thus  g^  (  j  ,  t,  a  )  is  the 
average  signal  received  on  the  i'th  hydrophone  at  time  t  from  the 
ensemble  of  target-medium  combinations  corresnonding  to  a  parameter 
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Fig.  3a  THE  IMPULSE  RESPONSE  FILTERS 


a.  3b 


THE  MODIFIED  IMPULSE  RESPONSE  FILTERS 


vector  a  in  response  to  a  unit  impulse  emitted  From  the  transmitter  at 
time  t  -  f  •  The  set  of  random  filters  described  by 


r(T.t,  3  )^r,(  T,tf  a)].  i  3  1 1  2...  A,  (Eq.  5) 


incorporate  the  effects  of  all  parameters  not  contained  in  3.  .  The 
output  functions  are  expressed  by  the  vector  valued  function 


3  (t,  a  ,  ?  )  =  [V  (t,  a .  rj  ], 


i  3  1 1  2...  A,  (c.q,5) 


where  we  have  written  r  for  r  ,  t,  a  J  and  r.  for  r.  [  <*  ,  t,  aj 


and  where 


si  (t,  a  .  r\)  =  s(t)  *  [gt  (  t  .  t»  a  )  +  ^  (  t  .  ti  a  fl  (Eq.  7) 

Figure  4  illustrates  the  shape  of  (  T  ,  t,  a  )  and  r\  (  T  ,  t,  a  ) 

as  functions  of  t  for  constant  t.  The  example  refers  to  a  linear 
receiving  array,  an  extended  target  and  a  multipath  transmission  medium. 


It  was  stated  in  the  introduction  that  decisions  about  the  target 

situation  were  to  be  based  upon  measurements  Dn  the  acoustic  field 

propagating  to  the  receiver.  In  the  following  we  shall  consider  a 

"measurement"  to  be  some  approximate  finite-dimensional  representation 

of  the  acoustic  pressure  signals  received  on  the  hydrophones  of  the 

receiving  array  during  one  ping  interval  T  i.e.  the  interval  between 

0 

two  transmissions.  We  shall  be  specific  in  our  choice  of  representation 
and  consider  a  measurement  to  be  the  amplitude  values  obtained  by 
sampling  the  acoustic  pressure  signals  received  on  the  hydrophones  of  the 
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EXAMPLE  CF  IMPULSE  RESPONSES 


receiving  array  at  a  rate  of  2B  samples  per  second  where  B  is  the 

bandwidth  of  the  sonar  signals.  He  assume  for  simplicity  that  the 

spectra  of  the  sonar  signals  have  been  shifted  down  into  a  lowpass 

band  of  width  3.  The  dimension  of  a  measurement  is  therefore 

K  =  2»A.3.T  and  the  dimension  of  an  echo  signal  vector  2«AtB*T 
o  9 

where  T  is  the  total  length  of  the  signal.  The  sample  values  of  a 
s 

measurement  arranged  in  a  suitable  sequence  form  the  elements  of  a 
measurement  vector  x  .  Applying  our  sampling  scheme  to  the  signals 

r  -4  ^  \  — ^  *  ,1  % 

s  (_t,  a  ,  r  J  of  equation  :3  we  obtain  the  random  signal  vector  3  [  a  J* 
For  a  given  g,  the  signal  vector  exhibits  random  fluctuations  from  one 
ping  intorva.  to  the  next.  The  fluctuations  are  caused  by  variations  in 

“4  “4  r  \ 

thG  random  rbsponss  functions  Gxprnssed  by  r  =  r  [  T  ,  t,  a 

In  later  chapters  is  will  ee  convenient  to  let  the  dimension  of  the 
signal  vector  be  equal  to  that  of  the  measurement  vectcr.  We  achieve 

this  by  assuming  the  signal  vector  to  have  one  element  for  eacn  sampling 

point  in  T  .  Elements  corresponding  to  samole  points  outside  the 
o 

signal  intervals  are  zero. 


3.  FLEf.'ENTG  CF  0ECI5IG’J  3V3TEV3 

The  tasiq  of  the  decision  unit  of  fig.  2  is  to  decide  which  targets, 
if  any,  are  present  in  the  surveyed  area.  Targets  are  identified  by 
their  parameter  vectors  a  .  Decisions  are  mads  at  the  end  of  each  • 
ping  interval  and  are  based  upon  the  measurement  x  received  in  the 

■4 

corresponding  ping  interval.  The  vector  x  contains  a  noise  component 

■4 

'n  comprising  all  ^n desired  signals  appearing  at  the  hydrophone  output 
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terminals  including  system  circuit  noise  referred  to  the  nydroohona 
terminals.  The  principal  components  of  n  will  be  ambient  noise, 
reverberations  and  in  some  cases  flaw  and  cavitation  noise.  If  a 
target  with  parameter  vector  a  is  present  a  corresponding  random 
signal  vector  s  (  a)  is  embedded  in  the  noise  n.  ’.Ve  shall  make 
the  simplifying  assumption  that  different  values  cf  a  are  mutually 
exclusive  i.e.  at  most  one  target  can  be  present  in  the  surveyed  area 
in  any  one  ping  interval.  This  assumption  does  net  alter  the  basic 
structure  of  the  signal  processing  device  of  the  sonar  receiver.  7/e 
shall  furthermore  assume  that  values  of  a  occurring  in  consecutive 
ping  intervals  are  statistically  independent.  This  is  not  strictly 
true  since  the  limited  speed  and  manouvring  capabilities  of  real 
targets  introduce  dependence  between  successive  a  values. 


Decision  processes  are  generally  cf  two  kinds:  detection  processes, 
which  decide  whether  any  signal  among  the  possible  signals  is 
present,  and  estimation  processes,  which  decide  which  signal  or 
signals  are  considered  to  be  present  when  a  detection  has  been  made. 
It  ma>  be  shown  that  two  types  of  inQut  information  are  required 
for  either  cf  the  two  orccesses.  Cne  is  the  conditional  orebability 
density  p  (a  /  x  )  expressing  the  probability  of  the  various 
target  vectors  a  after  the  measurement  x  has  been  acquired. 

For  the  detection  process  we  also  require  the  posterior  crobabiliey 
p  (  0  /  x  )  of  the  signal-absent  situation.  '.7a  obtain  this  from 
equation  1C  by  inserting  a  =  0.  The  other  input  to  the  decision 
device  is  a  cost  function  expressing  the  quality  of  decisions.  The 
cost  function  is  generated  by  the  user  of  the  system  and  may  vary 
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with  tin®.  In  this  aaoer  .vs  snail  cnly  bo  incorested  in  the  function 
p  [  g.  /  x  J  which  expresses  the  information  we  need  for  the  decision 

■4 

process  from  the  measurements  x  .  We  shall  proceed  to  study  this 

function  in  some  detail  and  shall  find  that  we  are  usually  only  able  to 

gen jrate  an  approximation  to  it.  Figure  5  illustrates  the  relations 

between  the  0-  ,  s  (  a  J ,  and  x  spaces.  The  transitions  from  th: 

-*  “♦  .  -4  . 

a-spaco  to  the  s  [  a  J-space  are  governed  by  the  probability 

density  function  p  s  (  a)/  a  J  which  reflects  the  statistics  of  the 

random  response  functions  described  by  r  =  r  (t  ,  t,  a  )  of 

equation  5,  The  transitions  from  the  5  [  a J-space  to  the  x-spaco  is 

r-+  . 

governed  by  the  probability  density  function  p  J^x  /  s  [  a  J  which 
depends  upon  the  statistics  of  the  noise  rt  .  uinco  the  noise  is 
additive  :n  can  write 


p|_x/s(a)J  =  PnL*~®^J  8) 

where  p  (.  n  J  is  the  probability  density  function  of  the  noise  n  . 

n 

The  probability  density  function  p  [  x  / a  j  governing  the  transition 

-4  -* 

from  the  a-space  to  the  :<  -space  can  now  be  written 


s 


where  the  integral  is  over  the  entire  s  (.  a  J-spacs.  Applying  the 
3ayes  Theorem  we  obtain  the  desired  function 


P  C 3/  ? )  -  p  till l1 

P  (  x  ] 


(Eq.  1C) 


Id 


Since 


shall  write 


P  ( 


is  a  constant  For  a  given  measurement  wo 


P  (  x  )  =  C_1 
and 

p  (  a  /  x  )  =  C.  p  (a),  p  (  ,<  /  a) 


(Eq.  11) 


/  ”4  »  x  4 

where  C  is  such  that  the  integral  of  p  (.  z/  x  j  w.r.t  a.  equals  one. 

r~*  \ 

The  function  p  [  Q.  j  is  the  probability  density  function  of  z  prior  to 
the  measurement  x  .  It  has  frequently  been  argued  that  the  function 
p  (  a)  is  difficult  ta  define  or  devoid  of  physical  meaning.  .7s  shall 
here  take  the  attitude  that  the  form  of  p  (a  ),  whether  of  the  minimax 
type,  constant  probability  type  or  otherwise,  must  be  chosen  by  the 
system  user  in  accordance  with  the  situation  he  believes  himself  tc  be 
in.  Presumably  the  system  user  is  the  best  source  of  prior  probabilities; 
he  is  in  any  case  the  only  one  available.  The  function  p  (  x  /  j) 
occurring  in  equation  1C  is  the  one  we  wish  to  obtain  from  the  processing 
unit  of  the  sonar  receiver  in  fig.  2.  This  function  will  be  the 
principal  theme  of  Chapters  %  and  5. 


77 e  remark  at  this  paint  that  equation  1 1  may  bo  extended  to  e  situation 

in  which  it  is  desired  to  make  decisions  based  cn  a  sequence  of  m 

-*  -»  -♦ 

consecutive  measurements  x  .  x„  ...  x  i.o.  to  -situations  in  which  it 

12  m 

is  desired  to  take  into  account  the  statistics;  dependence  between  z 
values  obtained  from  a  target  in  m  consecutive  sing  intervals. 
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4.  ADAPTIVE  PROCESSING 


The  structure  of  the  signal  processing  device  of  a  detection  system 
will  often  deoend  upon  one  or  more  parameters  characterizing 'the 
input  sic- a Is .  If  one  such  parameter  is  known  to  exhibit  random 
variations  with  time,  the  structure  of  the  processor  can  be  arrived  at 
in  two  ways.  If  the  variations  are  stationary  in  nature  and  their 
statistical  properties  are  known  wa  can  choose  that  processor  'which 
gives  the  best  average  performance.  Alternatively  we  may  carry  out 
measurements  in  order  to  estimate  the  true  value  of  the  parameter  and 
adjust  the  structure  of  the  processor  in  accordance  with  the  estimate. 
This  technique  reduces  the  uncertainty  about  the  unknown  parameter  and 
improves  the  system  performance.  The  latter  type  of  processor  is 
referred  to  as  an  adaptive  processor.  In  this  chapter  we  shall  study 
processors  which  adapt  to  a  randomly  varying  ncise  field. 


It  was  pointed  out  in  the  previous  section  that  the  information  required 

for  the  decision  process  from  the  measurements  "x  was  contained  in  the 
r  ^  ■n 

function  p  [  x  /  aj  of  equation  9.  The  discussion  about  adaptive 
systems  will  relate  to  the  structure  of  the  processing  device  from 
which  p  (  x  / gt  )  is  obtained.  V/s  reproduce  equation  9 


p  (  1  /  i]  -  j  P  p  /  3  c  3)  j  -D  [  ?  (3 )  /  3]  -d?  c  Si 


and  recall  that  the  function  p  |  "x  /  s  ( & )  j  under  the  integral  sign 
is  determined  by  the  probability  density  function  of  the  noise  vector 
n  in  accordance  with  equation  8,  and  that  p  |  "I  (a)  /a  J  reflects 
the  statistical  properties  of  the  targets  and  the  transmission  medium. 
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Thus  to  ub'tain  P  C  x  /  &  )  WB  require  a  noise-dependent  processing 
device  to  evaluate  p  £  ?  /  s  (a)J  for  each  3(a)  and  an  averaging 
device  depending  upon  p£s  (a)  /a  J.  These  are  shown  in  fig.  3.  Vie 
shall  proceed  to  study  the  ncise-depsndent  processor  and  shall  leave  eh 
discussion  of  the  averaging  device  to  section  5. 

The  discussion  will  proceed  under  the  assumption  that  the  noise  is 
gaussian.  This  implies  that  the  noise-dependent  processor  depends  only 
upon  the  second  order  statistical  properties  of  the  noise  field.  Vie 
remark  that  a  gauss.',  ^-equivalent  suboptimal  system  may  be  defined  from 
the  second  order  statistics  of  any  non-gaussian  noise  field.  Under  the 
gaussian  hypothesis  the  noise  vector  n  is  a  K-dimensionai  gaussian 
random  variable  with  probability  density  function  (Ref.  5  Chap.  24) 


p  (  n  )  ■  (2  Tt)  .  (Oet  wl)  .  exp  -  \  n  **  .  3  .  n 


(Eq.  ‘2) 


where  n  is  the  noise  column  vector  with  elements  n^,  i  =  1,  2...  K, 


n  is  the  transpose  of  n, 

■4  ->_'l 

3  =  f.I  is  the  inverse  moment  matrix, 


M  =  j  m.  .1 
L  ljJ 


"l  _*  r 

is  the  moment  matrix  of  n  with  m.  .  =  :  |  n,<n,  , 

ljJ  1J  L  1  J  J 


Det  M  -  the  determinant  of  M 

We  note  at  this  point  that  the  matrices  >V1  and  3  are  real  and 

*4  4  t  -*  4t 

symmetrical  such  that  M  3  M  and  3=3  .  3y  equation  8  we  are 

now  able  to  write 


IB 


Pp/sC3?)]  POT/ST) 


i 


_ 

p[x  /  t  C  a)  ]-  (2  r)  r.  (Det  Ml-4  exp  -  4  {f  -  t*)  .3.  (x  -  3) 


D  .  exp  (x  .  J  .  3  -  i  5"  .  3  .  s)  (Eq.  13) 
n 


where 


D  = 


_ 

(2TT  )  T  .  (Det  M)“  5  ^  .  Q  .  x 


and  where  for  simplicity  we  have  written  "s  for  s  (  a)*  Fig.  7 
illustrates  the  formation  of  the  exponent 


Z  “*  X  *  di  1  3  "  V  S  •  uf  •  o 


(Eq.  U) 


of  equation  13.  In  many  oractical  situations  the  elements  of  the  matrix 
3  =  I  q.  .  will  exhibit  random  variations  with  time.  A  system  which 

-  U  J 

keens  track  of  the  elements  q .  .  and  adjusts  the  processors  F  and 

^  J  * 

of  Fig.  7  in  accordance  with  their  values  is  the  adaptive  system 
we  are  seeking.  For  reasons  which  will  become  clear  later  we  shall 
refer  to  this  as  a  regression  type  system.  3efore  we  study  this  in 
detail  similar  adaptive  systems  based  uocn  second  order  noise 
statistics  will  be  reviewed  briefly. 


Three  types  of  signal  filters  applicable  to  detection  and  extraction 
systems  have  received  much  attention  in  recent  years.  They  are 

(1)  the  maximum  signal-to-noise  ratio  filter  (Ref.  11  ) 

(2)  the  minimum  signal  distortion  or  Wiener  filter  (Ref.  3  )  and 

(3)  the  regression  type  filter  (Ref.  2  ). 
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Whan  ussd  to  detect  monochromatic  signals  the  spatial  filtering 
properties  of  these  filters  have  been  shown  to  be  identical  (Ref.  S  ). 
Their  temporal  filtering  properties  are,  however,  different,  as  are  also 
their  areas  of  application.  In  principles  the  three  filter  types  can 
operate  in  an  adaptive  mode.  An  adaptive  "maximum  signal-to-noise  ratio" 
filter  has  bean  described  by  Shor  (Ref.  12  )  and  an  adaptive  Wiener  filter 
by  Widrow  at  al  (Ref.  13  ).  This  last  article  is  the  principal  reference 
of  this  chapter.  The  adaptive  technique  suggested  for  the  regression 
filter  to  be  presented  is  taken  directly  frcn  oho  wore  of  ,/idrow  at  al. 

In  order  to  expose  the  salient  features  of  the  adaptive  regression 

system  and  avoid  undue  compexity  wa  shall  assume  that  the  sonar  system  is 

ambient  noise  limited  and  that  the  variations  in  q.  .  are  small  over  c.ns 

ij 

interval  T  .  The  reverberation  field  associated  with  a  moving  sonar 

0 

platform  is  considered  to  exhibit  too  rapid  variations  to  bo  amenable  to 
adaptive  techniques.  For  a  stationary  platform  the  reverberation  field, 
although  rapidly  varying  within  each  sing  interval,  may  exhibit  sIdw 
variations  from  ana  ping  intarval  to  tha  next  when  values  at  identical 
positions  within  tho  ping  intervals  are  considered.  Adaptive  techniques 
may  therefore  be  applied. 


Proceeding  with  the  ambient-noisa-limitod  case  we  turn  to  fig.  0  which 

illustrates  a  composite  measurement  y  consisting  of  the  j  measurement 

m  * 

-4  -4 

vectors  preceding  ths  vector  x  •  Thus  y  contains  the  input  data 

m  m 

received  in  an  interval  of  length  T  3  j.T  ,  We  shall  think  of  T 

y  °  y 

as  an  interval  over  which  3  remains  essentially  constant.  Provided 

-4 

no  signals  are  present  in  ths  measurement  y  and  provided  the 

■4  m 

variations  in  ths  elements  of  the  matrix  „  are  stationary  with  :nc.vn 
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( 


statistical  properties  we  can  obtain  a  posterior  orccability  density 

-4  -4  -* 

function  p  (Q  /  y  )  of  U.  In  principle  this  should  b ;  jeed  to  obtain 
m 

the  conditional  forward  probability  function 

-  r V 

P  !~x  /  s  (  a),  y  P  p/s  (a),  jj  .  p  (J  /  y]  .  dd  (Eq.  15j 


which  takes  the  place  of  p  j^x  /  s*  (  a)  J  in  equation  9.  In  many  cases 
the  statistical  properties  of  Q  are  unknown  such  that  P  [  Q  /  y  J 
can  not  be  obtained.  Vie  can  then  choose  to  make  estimates  of  J  based 
on  the  composite  measurement  vector  y  and  use  the  estimates  in 
equation  11  to  generate  z.  The  corresponding  processor  of  Fig.  7  is 

-4 

then  adapting  to  a  changing  matrix  Q.  Vie  shall  study  two  techniques  fer 

*4  4 

estimating  tha  matrix  J  from  the  measurements  y.  strictly  speaking 

we  require  that  the  y’  s  do  not  contain  any  signal  components.  Thus, 

prior  to  each  decision  interval  (corresponding  to  x  in  Fig.  6),  we 

m 

require  a  learning  phase  of  length  T  in  which  the  characteristics  of 

y 

the  noise  field  are  assessed.  This  feature  appears  to  be  common  to  all 

noise  adapting  systems.  In  real  systems  the  decision  interval  is  pushed 

forward  along  the  time  axis  an  amount  T  at  the  beginning  cf  each  ping 

o 

interval.  The  noise  adaption  process  will  therefore  proceed  smoothly  un 
to  the  point  when  a  target  signal  appears  in  the  received  waveforms.  From 
then  and  onwards  the  adaptive  processor  will  consider  the  target  sign il 
as  part  of  the  noise  field.  The  extent  to  which  it  will  succeed  in 
suppressing  the  signal  depends  upon  the  angular  velocity  of  the  target 
relative  to  the  settling  time  of  the  prccessor.  (The  settling  time  is 
the  tine  required  for  the  adaptive  prccessor  to  resnend  tc  a  change  in 


the  -ambient  noise  field.] 


The  first  tyoe  of  adaptive  technique  will  be  touched  upon  very 

lightly.  It  is  the  direct  aporoach  of  estimating  J  through 

■4 

estimates  of  the  moment  matrix  M.  Thus  estimates  of  the  elements 

-4  -4 

m,  .  of  M  are  obtained  in  each  learning  interval  T  whereupon  Q 

ij  y 

is  obtained  through  matrix  inversion.  If  the  power  spectra  of  the 
received  noise  waveforms  and  the  power  spectra  of  the  variations  in  the 

moments  m.  ,  sire  known  it  is  a  simple  task  to  design  good  estimators 

•  J 

for  the  moments  m  .  This  technique  is  rather  impractical  since  it 

requires  a  large  number  of  correlators  for  determination  of  the  moments 

m,  .  and  since  it  involves  the  operation  of  matrix  inversion, 
ij 


The  second  type  of  adaptive  technique  is  based  upon  the  gradient-search 
method  discussed  in  Ref.  13  .  In  order  to  understand  Phis  method  it  is 

-4 

necessary  to  paint  out  certain  properties  of  the  matrix  ^  and  the 
vector 


v  *  Q  .  x 


Vh 


K 

V 

/ 

i  =1 


q,  .  .  x 

:u  i 


with  elements 


h  =  1,  2,  ...  K 


-6) 


appearing  in  the  expression  for  z  in  equation  14.  Consider  first  the 
case  of  a  time-invariant  gaussian  noise  field  with  an  associated 
K  x  K  dimensional  symmetrical  noise  moment  matrix  M  with  inverse  Q, 

-4  *4 

Assuming  noise  only  to  be  present,  i.e.  x  =  n,  it  may  be  shown  that 
v.  of  equation  15  is  equal  to  the  residue  rr.  of  the  h'th 


hi 


observation  element  x  in  x  divided  by  the  variance  c  ri 

h  h 

this  residue  (Ref.  5  Chap.  23).  Thus 


2  >  of 


S  7  << 


1, 


(=R.  17) 
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1 


I 


It  is  important  tc  grasp  tha  physical  significance  of  the  residue  ri  . 

h 

It  represents  the  difference  between  the  observation  element  x  and 

h 


the  best  linear  estimate  x  of 

h 

K-l  elements  of  x.  Thus 


obtained  from  the  remaining 


n  =  (x  -  x  ) 
h  h  h 


(Eq.  IB) 


and  where  the  constants  3,  are  adjusted  to  give  minimum  of  the 

residue,  variance  <  Tl2  >•  (The  symbol  T,  f  h  under  the  summation 

h 

sign  indicates  summation  ever  all  £  oxcart  i  =  h.)  By  the  gradient 

search  technique  the  regression  coefficients  q  are  adjusted  in 

small  steps  in  a  systematic  way  such  that  the  minimum  value  of  <  r2  •> 

h 

is  approached  along  the  path  of  steepest  descent.  To  obtain  the 

desired  quantities  v  of  equation  17  the  residues  n  ,  thus 

h  1  h 

obtained  must  be  divided  by  estimates  of  the  residue  variance  <r  r.E  > 

h 

The  gradient-search  method  will  now  be  explained  in  some  detail. 


Referring  to  equation  13  we  shall  represent  the  set  of  coefficients 

S  ,  associated  with  the  estimate  x  by  the  vector  £ .  .  Thus 
hT  h  n 

,  l  =  1 ,  2,  •••  (h-i),  (h+1 ;  ...K  (Eq .  19) 


Since  will  change  with  time  we  shall  write  p"^  (j)  for  the 
value  of  £  in  the  j'th  ping  interval.  By  the  gradient-searc.1- 
technique  we  apply  zhs  followine  recurrence  formula  to  3. 
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9h  Cj  +  i)  -eh  tj)  <■  k  •  7 


\  ^ 


[Eq.  -C ' 


where 


7  <r 


2 

h 


k 

s 


(j) 


=  negative  scalar  constant  controlling  the  rate  of  convergenc 
and  stability. 

>=  gradient  vector  of  residue  variance  with  respect  to  3  • 


In  practical  situations  we  do  not  have  access  to  the  true  gradient 

vectors.  Following  Ref.  13  we  shall  use  the  gradient  of  ri2  (j) 

1  h 

as  estimates  for  7  <  ri2  (j)  >•  From  equation  18  this  gradient  is 

h 


7  rfh(j)  -  -  2  rh(j).  3h(j)  fa.  20 

where  x"*(j)  is  the  observation  vector  corresponding  to  the  j’th 
h 

ping  interval  with  the  h'th  element  deleted.  Substituting  this 
result  for  the  gradient  vector  cf  equation  20  we  obtain  the  final 
recurrence  formula 


a.  (j+o  -  B.  (j)  -2k 

h  h  s 


Cj) 


’h(j) 


(Eq.  22) 


Widrow  et  al  show  that  the  estimates  7  ri2  f  \  )  are  unbiased 

1  h 

and  that  the  elements  of  3  converge  in  the  mean  towards  cha 

h 

true  regression  coefficients.  Figure  9  illustrates  the  adaptive 

formation  of  the  quantities  v  of  equation  17.  The  estimate  of  the 

h 

residue  variance  <T12  fj)>  is  a  suitably  weighted  sum  of  the  square  cf 

h 

past  values  of  ri .  •  The  number  of  terms  in  the  summation  depends  on 
h 

the  rate  of  variation  cf  the  noise  field. 
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.Updating  of  :  Ai  ( j*l )  =  ^h(j )- 2* ks *??h(j )*Xh(j ) 


THE  ADAPTIVE  PROCESSOR 


Up  to  this  point  the  discussion  of  the  adaptive  regression  system  has 

bean  cantrad  upon  the  residue  of  a  single  element  of  x  .  At  a  first 

glance  it  may  appear  that  a  seoarata  adaption  process  must  be 

instrumented  for  each  of  the  K  residues  of  the  observation 

vector  x.  Vie  shall  show,  howavar,  that  one  adaptive  processor  for 

each  hydrophone  channel  is  sufficient.  To  see  this  we  change  our 

point  of  view  slightly  and  consider  the  regression  processor  as  a 

digital  filter  applied  to  the  samplad  input  waveforms.  Fig.  10 

illustrates  the  samples  belonging  to  an  observation  vector  x 

arranged  in  A  discrete  time  series,  one  for  each  hydrophone.  Vie 

focus  attention  on  the  observation  alement  x,  which  occupies  the 

h 

central  position  among  the  samples  from  the  i'th  hydrophone.  Next 

we  define  an  adaptive  digital  filter  with  weights  3  by  the  relation 

n</  . 

K 

=  I  3h-t*  \  ^Eq* 2 ^ 

i=  1 


where 


■p  =  residue  of  observation  element 

h 


=  - 1 


it  h,  are  the  variable  regression  coefficients  which  a re 


hh 

8  hi  ’  ^  ■ 

the  elements  of  g^  of  equation  19  . 


Equation  23  is  obtained  from  equation  16  by  introducing  =  -  ;  . 

The  filter  defined  by  equation  23  is  the  adaptive  regression  filter 
for  the  waveform  from  the  i'th  hydrophone.  As  new  samples  are 
shifted  into  the  filter  successive  samples  of  the  residue  waveform 
from  the  i'th  hydrophone  are  generated  at  the  output. 
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X 


Xh 


/ 

%=  /M-x, 

K 

=  ’  Xl  5  ^hh  ~  ~  ^ 


Fig.  10  ILLUSTRATION  TO  THE  FORMATION  OF  ^  FROM  x 


30 


In  changing  from  the  purely  vector-oriented  regression  processor  of 

equation  18  to  the  Filtor-oriented  processor  associated  with 

equation  23  we  have  implicitly  changed  our  concepts  about  the 

observations  x.  Originally  the  observations  were  generated  at  the 

rate  of  one  per  T  seconds,  i.e.  one  per  Ding  interval.  For  the 
a 

filter-oriented  processor  the  observations  change  by  A  new  samples 
introduced  and  A  old  samples  discarded  in  each  ping  interval.  In 
order  to  bring  out  clearly  the  filter  point  of  view  we  modify  the 
terminology  of  equation  23  as  follows 


l  -1 


ht 


i  -  1,  2,  ...  K 


(Eq.  24) 


where 

T|^  is  the  output  of  the  i*th  filter  in  the  h'th 

h 

sampling  interval, 

x  is  the  observation  element  occupying  position  "  l”  in 

nt 

the  observation  centred  on  sampling  time  "h",  and 
0  ^  is  the  filter  weight  occupying  position  "l"  among 
the  filter  weights  of  the  i'th  filter. 


At  this  point  we  mention  that  in  real  systems  we  shall  expect  that 


only  a  limited  number  of  observation  elements  x  in  the  neighbourhood 

hi 

of  sampling  time  "h"  will,  have  significant  influence  upon  y 

The  time  interval  embraced  by  the  digital  filter  response  may  then  be 

considerably  smaller  than  T  and  the  number  of  filter  weights 

o 

considerably  smaller  than  K. 
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A  bank  of  regression  filters,  one  for  each  hydrophone  and  each  filter 
followed  by  a  variable  gain  device  which  affectively  divides  the 
filter  output  by  an  estimate  of  the  output  variance,  make  up  the 
complete  adaptive  regression  processor.  These  are  shown  in  Fig.  11. 

The  structure  of  each  processor  is  as  shown  in  Fig.  9.  7/e  note  that 

a  second  set  of  identical  filters  are  required  for  the  signal-dependent 

-*  t  -*  -* 

tern  s  .  u  .  s  of  equation  14. 


5.  THE  SIGNAL  DEPENDENT  PnCCZGGCP 


We  recall  that  the  information  required  for  the  decision  process 

,  -4  \ 

from  each  measurement  x  was  contained  in  the  function  p^x  /  a  J 
where  a  can  vary  over  the  set  of  possible  target  parameters. 
From  equation  9 


p(  *  /  a] 


s 

j  o  [?/1  (2)  ] .  p  [2  (3)  /J 


In  the  previous  section  we  studied  the  adaptive  noise  dependent 
processor  of  a  noise  limited  sonar  receiver  operating  in  a  gaussian 
ambient  noise  field.  This  processor  evaluated  the  function 
(equations  13  E-  17) 


.  s  [  a  . 
(Eq.  25) 


appearing  under  the  integral 
s  (  a) .  In  this  section  we 


signal  Gf  equation  5,  as  a  function  of 
shall  study  the  structure  of  the  signal 
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dependent  processor  which  evaluates  P  (  x  /  oD  when 
p£  x  /  s  (  a)  is  as  given  above  and  p£s(a)/aj  i-3  known. 

Generally  speaking  the  integral  expression  for  p  (  x  /  ^ ) 
given  above  cannot  bo  developed  further  except  when  p^3  (a)  /aj 
has  a  simple  mathematical  form,  c.g.  a  multivariate  gaussian 
probability*  density  function,  cr  -when  the  exponent 


-*t  ,  -» 

z  »  s  (a 


(?)  .  3  .  ?*  (a)  •  3  .  3  (  S) 


in  the  exarassion  for  p^x/afa]!  lb  almost  always  small 
ccmparcd  to  unity  such  that  the  first  terms  of  a  power  series 
expansion  may  bo  used  to  represent  p  p  /  s  (  q.)  J  •  V/o  shall 
study  the  structure  of  the  processors  obtaineo  in  chose  special 
cases  and  comment  upon  the  syeo  of  signal  stacistics  required  for 
their  roa.icaticn. 

We  consider  first  the  case  when  s  ( ^ }  ic  composed  of  a 

•4  0  -4  ^ 

non-random  port  a  v  a  ;  and  an  additive  gaussian  perturbation 

o 

,  “♦  f  ~*y 

vector  u  [  a,1-  inus 

slai*3  l»j+uka( 

^  *4 

The  perturbation  vector  is  associated  with  a  moment  matrix  U  [  a' . 
(As  mentioned  in  Chapter  2  the  number  of  elements  In  the  signal 

*  -4  % 

vector  3  l  aj  is  considered  to  be  equal  to  the  number  of  elements 

•4 

in  x  .  The  matrix  U  [  aj  is  therefore  also  a  K  x  K  symmetrical 
matrix. }  3inoo  u  (  a}  is  gaussian  wo  can  write  down 


J 


immediately 


-1 


p(x  /  ft  )  =  □(  g]  .  exp  -  i[  X  fc  -  sofc  (  2)  !  .  [m  +  U  C  a]  j  .  p  -  SQ  (3) 

(Eq.  25) 


with  0  (a)  =  [2tt  )" 


K 
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Dat  Cm  +  U  (  a))]  "  ^ 


W 8  can  split  equation  25  up  into  a  product  of  the  following  two 
components: 

(l)  the  U  (a)  dependent  term 


=  D  (  a  )  exp  -  i  2fc  L«  +  U  (o) 


-  1 


and 


the 


u  (3) 


and 


•0(SJ 


dependent  term 


<=2  =  exp 


’♦t 

v 


-•t  /•-*  i 

-  sq  (a)  . 


I  + 


,7  ,  -*> 
u  l  a; 


If  the  energy  in  the  signal  perturbation  component  tends  to  zero,  the 
first  term  will  tend  to  a  constant  whereas  the  second  tend  to 
exp  z  with  (see  Fig.  7  and  equation  1  l) 


z 


We  see  that  z  contains  the  correlation  term  between  the  observe sion 

as  it  appears  at  the  output  of  the  noise  dependent  processor  and  the 

-*  ,  -»  , 

average  signal  3  [  g.  J .  If,  however,  the  mean  signal  tends  to  zero 

o 

we  must  work  with  the  first  term.  When  the  elements  of  U  (  a)  are 
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small  compared  to  the  slemants  of  M  i.o.  for  small  incut  signal-to- 
noise  ratios,  we  can  writa 

C  =>  D  exp  -r.  v  .  U  (  a)  .  u 
1  n 

with 

D  =  (2r  )  3  (Oet  !;)  exp  -  x***  u  x* 

n 

The  exponent  v  .  U  [  aj  .  v  can  be  shown  to  be  equivalent  to  the 

output  of  an  optimum  energy  detector  searching  for  noise-like 

signals  from  a  target  with  parameters  a.  This  situation  reoresents 

one  of  extrema  signal  distortion.  As  would  be  expected  from  the 

assumption  of  a  gaussian  signal  perturbation  vector  equation  2  shows 

that  the  statistical  knowledge  required  about  the  signals  is  the 

average  value  3  [  a  J  and  the  second  order  moments  between  the 

o 

elements  of  the  signal  perturbation  vector  u  (  aj.  These  quantities 
must  be  known  as  functions  of  the  target  parameter  a  • 


The  second  special  case  to  be  studied  relates  tc  the  situation  where 
the  standard  deviation  of 


z 


is  much  smaller  than  unity  such  that  z  itself  is  almost  always, 
small  compared  to  unity.  'Us  can  then  use  the  first  three  terms  cf  a 
power  series  expansion  to  represent  pl”x/"s(g)  i  of  equation  '.2 
and  obtain  when  terms  involving  third  and  higher  order  signal  moments 
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9ts  discarded 


p(* /a)  «on  +  °n  .  a  .  s*  (  a)  j 

t  on  [v c .  u  (  3)  .  ?  +  i  .  to  (3  ))s ]  (Eg.  37) 

Since  we  have  discarded  all  terms  involving  signal  moments  higher 
than  the  second,  the  processor  is  again  described  in  terms  of  first 
and  second  order  signal  moments. 

A  large  amount  of  present-day  research  on  the  acoustic  propagation 
properties  of  th8  ocean  relates  to  the  study  of  the  average  signals 
and  the  second  order  time-space  moments  of  the  perturbation  components 
of  signals  transmitted  over  direct,  surface-reflected  and 
bottom-reflected  paths.  The  relations  between  these  signal  components 
and  the  physical  properties  of  the  ocean  and  its  boundaries  are  of 
primary  importance  in  such  studies.  It  should  be  noted  that  such 
studies,  although  already  vary  complex,  are  not  complete  since  the 
real  sonar  situation  involves  the  transmission  of  signal  energy 
to  and  from  a  complex  extended  reflector.  The  information  required 

•  f  “4  \ 

about  the  medium  and  targets  to  obtain  the  quantities  s  (.  a  J  and 

o 

r~*\ 

U  (a  J  on  which  the  signal-dependent  processor  depends  can  be 
related  to  the  target-medium  impulse  responses  introduced  in 
Chapter  2.  Speaking  in  terms  of  analogue  waveforms  rather  than 
sampled  waveforms  the  average  signal  received  on  the  i'th 

t  -4 

hydrophone  from  a  target  with  parameter  a  in  response  to  a 
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transmitted  signal  s(t)  is 


00 


t,  ft 


>} 


(t  -  T  ) 


d  t 


(Eq .  2e  ' 


Ths  average  signal  is  seen  tc  depend  upon  the  function 


introduced  in  chapter  2  and  referred  to  in  Fig.  3b  and  i.  A 

■4  .  *4 « 

typical  element  of  the  signal  perturbation  matrix  U  (.  g.J  will  be 
the  average  cross-product  between  the  oerturbaticn  signals  received 
an  hydrophones  i  and  j  at  times  t^  and  t^.  (Note  that  t^  ? 
t  specify  positions  within  a  ping  interval.)  Thus  writing 
u  C i f  j»  t ^ ,  t  )  for  the  matrix  element 


00  OS 


0  V 

o  o 


t2;  a) 


S  (t2  -  t2)  .  d-r  d,2 
s(t1  -  Tl )  •  -  t2)  .  d 
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1  ? 


(Eq.  32) 
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wh  era 


^  [t,  a,  r\^=  e|V  ri)  T  and 

Rij  (T1  •  t2;  t1 ,  t2;  3  )  =  E  J"  r.  (Tif  t1 ,  a  )  .  r.  (r  ^  a  )j- 

We  see  that  the  matrix  elements  depend  upon  the  correlation  function 
^ij(Tl'  T?’  S’  fc2*  ^  between  the  random  components  of  the  impulse 
responses.  In  Chapter  2  wo  assumed  the  target-medium  impulse  responses 
to  be  stationary  random  functions  of  time  t.  The  correlation  function 
will  then  only  depend  upon  the  time  difference  t  -  t^.  Finally,  if 

*4  -♦ 

p  and  p  are  the  position  vectors  of  the  i'th  and  i'th 

i  .] 

hydrophones  relative  to  the  receiver  reference  point  we  can  write  the 
time-space  correlation  function  between  two  impulse  responses 


P . 
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5.  C0NCLL13ICM3 


The  task  of  an  active  sonar  system  is  to  survey  an  ocean  area  in  the  sense 
of  reporting  upon  the  movements  of  targets  entering  this  area.  The 
output  of  the  system  is  decisions  about  what  is  considered  to  be  the  true 
target  situation  in  the  surveyed  area.  The  decisions  are  based  upon 
measurements  obtained  from  the  acoustic  signals  received  on  the  hydrophones 
of  the  receiving  array.  We  define  a  "measurement"  to  be  a  suitable  finite 
dimensional  representation  of  the  acoustic  waveforms  received  in  one  ping 
interval.  Apart  from  the  acoustic  transmission  properties  of  the  ocean 
and  the  reflecting  properties  of  the  targets,  the  characteristics  of  the 
measurements  depend  upon  the  transmitted  signal,  the  array  geometry  and 
the  signal  representation  chosen,  all  of  which  are  under  the  control  of 
the  system  design :r. 

'.Vo  consider  the  optimum  structure  of  a  sonar  receiver  from  the  hydrophones 
to  the  input  terminal  of  the  decision-making  device  for  a  given  measuring 
scheme  i.o.  for  a  given  transmitted  signal,  array  geometry  and  signal 
representation.  We  note  that  the  complementary  problem  of  deciding  uocn 
the  measuring  scheme  which  will  give  optimum  system  performance  for  a 
given  set  of  performance  criteria  is  far  more  complex. 

There  are  many  parameters  which  influence  the  echo  signals  received  from 
targets  in  response  to  a  transmitted  signal.  3ome  of  those  relate  to  the 
movements  and  structures  of  targets  and  some  to  the  acoustic  orepagstion 
properties  of  the  ocean  and  the  reflecting  properties  of  its  boundaries. 

A  sonar  system  is  designed  to  provide  information  on  seme  of  these 
parameters,  usually  these  directly  connected  with  the  decisions.  The 
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remaining  properties  ar  ted  as  random  parameters  from  the  point  of 

view  of  the  sonar  receiver,  which  must  be  constructed  so  as  to  be 
optimum  in  the  average  sense  against  these  unknown  parameters.  Generally 
speaking  system  performance  and  complexity  goes  up  with  the  number  of 
parameters  included  among  the  search  parameters  of  a  sonar  receiver. 

Target  echo  signals  arc  always  received  in  the  presence  of  interfering 
noises.  The  most  important  sources  of  noise  are  ambient  noise, 
reverberations  and  flow  noise.  The  interfering  noise  field,  which  at 
best  may  be  stationary  and  unknown,  will  generally  exhibit  variations  with 
time.  In  situations  involving  varying  and  strongly  non-isotropic  ncise 
fields  considerable  improvements  in  processing  gain  and  system  performance 
may  be  obtained  by  the  application  cf  adaptive  processing  techniques.  Ve 
consider  the  application  of  such  techniques  to  slowly  varying  gaussian 
ambient  noise  fields.  They  are  in  principle  also  applicable  to  the 
reverberation  field  from  a  stationary  sonar  platform,  although  with 
considerable  increase  in  system  complexity,  and  to  flow  and  cavitation 
noise  fields  which  remain  stationary  over  time  intervals  long  enough  to. 
permit  adaptation.  The  adaptive  processor  for  gaussian  noise  depends 
only  upon  second  order  noise  moments.  Gaussian-equivalent  subcctimum 
adaptive  systems  may  be  designed  upon  the  second  order  moments  of  ary 
non-gaussian  noise  field. 

It  was  mentioned  above  that  the  family  of  echo  signals  appearing  at  the 
hydrophones  of  a  sonar  system  depend  on  a  large  number  of  parameters 
relating  to  the  characteristics  of  the  targets  and  tho  ocean.  The  echo 
signals  are  described  in  terms  of  some  of  these  parameters,  the 
"search  parameters",  the  remaining  ones  being  treated  as  random  parameters. 
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Xhus  for  each  sat  of  values  of  the  "search  parameters"  there  exists  an 
ensemble  of  possible  echo  signals  and  the  corresponding  receiver  channel 
must  be  optimum  in  the  average  sense  against  this  signal  ensemble.  A 
practical  useful  formulation  for  the  average  receiver  channel  can  only  be 
obtained  for  some  simple  cases,  e.g.  when  the  echo  signals  consist  of 
gaussian  amplitude  perturbations  on  some  mean  signal  or  when  the  input 
signal-to-noise  ratio  is  low  enough  to  permit  a  description  of  the 
channel  in  terms  of  first  and  second  order  signal  moments.  The 
description  of  the  receiver  channel  for  these  cases  contains  terms 
indicating  a  cross-correlation  between  the  received  waveforms  and  the 
mean  signal,  and  energy  detection  of  the  signal  perturbation  comnonant. 
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